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Abstract 

This is a pedagogical account on reaction-diffusion systems and their 
relationship with integrable quantum spin chains. 

Reaction-diffusion systems are paradigmatic examples of non-equi- 
librium systems. Their long-time behaviour is strongly influenced 
through fluctuation effects in low dimensions which renders the ha- 
bitual mean-field kinetic equations inapplicable. Starting from the 
master equation rewritten as a Schrodinger equation with imaginary 
time, the associated quantum Hamiltonian of certain one-dimensional 
reaction-diffusion models is closely related to integrable magnetic chains. 
The relationship with the Hecke algebra and its quotients allows to 
identify integrable reaction-diffusion models and, through the Baxter- 
ization procedure, relate them to the solutions of Yang-Baxter equa- 
tions which can be solved via the Bethe ansatz. Methods such as 
spectral and partial integrability, free fermions, similarity transforma- 
tions or diffusion algebras are reviewed, with several concrete examples 
treated explicitly. An outlook on how the recently-introduced concept 
of local scale invariance might become useful in the description of non- 
equilibrium ageing phenomena is presented, with particular emphasis 
on the kinetic Ising model with Glauber dynamics. 



Ho be published in: A. Kundu (ed) Classical and Quantum Nonlinear Integrable Sys- 
tems: Theory and Applications, Institute of Physics Publishing (Bristol 2003) 
ISBN 07503 09598 



1 Reaction-diffusion processes 



The understanding of non-equilibrium statistical physics is still much more 
incomplete than that of equilibrium theory, due to the absence of an analogue 
of the Boltzmann-Gibbs approach and in spite of considerable recent progress 
PP. Therefore, non-equilibrium systems have to be specified by some defin- 
ing dynamical rules which are then analyzed. The topic has received a lot of 
attention and many reviews exist, e.g. [2HHJIUI2JIIEII1IE1III]- Exactly solvable 
systems far from equilibrium have been recently reviewed in a nice way ^H] ■ 
Here I present a pedagogically-minded introduction to the application of a 
few standard tools from one-dimensional integrable quantum systems to non- 
equilibrium statistical mechanics. After recalling why standard descriptions 
such a kinetic or reaction-diffusion differential equations are in general insuf- 
ficient in one dimension, we remind the reader in section 2 on the quantum 
Hamiltonian formulation of non-equilibrium processes which in turn is based 
on the master equation. Section 3 recalls a few basic facts about Hecke 
algebras. These building blocks are used in sections 4 and 5 to show ex- 
plicitly the integrability of certain single-species reaction-diffusion processes, 
through their relation to integrable vertex models. Section 6 reviews some 
further methods such as spectral and partial integrability, the free-fermion 
technique, similarity transformations or diffusion algebras. We close in sec- 
tion 7 with an outlook on how the recently-introduced concept of local scale 
invariance might become useful in the description of non-equilibrium ageing 
phenomena. 

I have made no effort to provide a complete bibliography. This may be 
found in the excellent reviews quoted above. 
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Figure 1: Microscopic reactions in the diffusive pair-annihilation process. 

A class of non-equilibrium models which are particularly simple to for- 
mulate are the so-called reaction- diffusion processes. Consider the following 
example: particles of a single kind (species) move on a lattice (figure^). Each 
site of the lattice may be either empty (denoted by o) or else be occupied 
by a single particle (denoted by •). The particles are allowed to undergo the 
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following movements, see figure ^ which involve the states of two nearest- 
neighbour sites: 



The first of the allowed movements in (0) is reversible while the second is 
not. A typical question is then for the long-time behaviour of quantities such 
as the mean particle density n(t). Trivially n(t) decreases with increasing 
time t but different long-time asymptotic behaviours such as n(t) ~ t~ a or 
n(t) ~ e~ l l T are conceivable. The oldest approach to this problem was intro- 
duced by Smoluchowski [TT] and consists of writing down kinetic equations, 
e.g. for the spatially averaged density n(t) and one obtains d t n(t) = —\n(t) 2 
for the problem at hand, where A = 4a. With the initial condition n(0) = n , 
the solution n{t) = rto(l + ra At) _1 ~ (At) -1 is easily found and apparently 
answers the physical question. A slightly more involved version of this argu- 
ment does allow for spatial variation of the density n = n(r, t) and considers 
a reaction- diffusion equation 



While the analysis of such non-linear partial differential equations is a for- 
midable problem in its own right, these equations do not yet capture the 
essential physics in low-spatial dimensions, as we now show. Rather, they 
must be considered as an approximation of mean-field type. 

In order to understand the approximative nature of equations such as (J2J), 
and following ^2], consider the mean particle-density in a large volume V 



o +• <-» • + o diffusion, with rate D 

• + • — > o + o pair-annihilation, with rate 2a 



(1) 




(2) 




(3) 



It then follows 




(4) 
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In the second line, GauB' theorem was used where dcr is normal to the surface 
dV (the flow through dV vanishes for large volumes V — > oo). The last 
line follows from the Cauchy-Schwarz inequality. Together with the initial 
condition n(0) = no, the inequality d t n{t) < —Xn(t) 2 yields the bound 

n{t) < U ° < 1 (5) 
w ~ 1 + n Xt ~ Xt K J 

for all times t > 0. However, the model defined above can be solved exactly 
in one spatial dimension (in a setting defined precisely in section 3), provided 
D = 2a. The exact mean particle-density is given by 



n(t) = n e 

1 



ADt 



I (4Dt) + 2(1 - n ) - 2n ) k - 1 I k (4Dt) 



k=l 



V8nD 



r 1/2 ; t -> oo (6) 



where Ik is the modified Bessel function of order k. Clearly, for large times 
the exact mean density n(t) decreases considerably slower than even the 
upper bound (J3J) derived from the reaction-diffusion equation (J2J. 

The failure of eq. (0) to describe correctly the long-time behaviour can be 
understood from the following heuristic argument In the long-time limit, 
the particle-density should already be very low and it is conceivable that at 
most one annihilation reaction takes place at a given time. Let L = L(t) be 
the typical distance between two particles. Then the time needed by diffusive 
motion to overcome this distance is of the order t ~ L{t) 2 - On the other hand, 
the mean particle density is n(t) ~ L{t)~ d in d spatial dimensions and this 
argument would give n(t) ~ t~ d l 2 . Therefore, the assumption implicit in 
equations such as (J2J that diffusive motion can render the system sufficiently 
homogeneous fails in low dimensions (in our model for d < 2) and one rather 
has the long-time behaviour 



Therefore d* = 2 is the upper critical dimension of the diffusive pair-annihi- 
lation process. For dimensions d > d*, reaction-diffusion equations should be 
expected to give qualitatively correct results and entire branches of physical 
chemistry are built on this. On the other hand, for low-dimensional structures 
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Table 1: Measured decay exponent y of the mean exciton density n(t) ~ t~ y 
on polymer chains. The error bar for TMMC comes from averaging over the 
results of |T7j for different initial particle densities. 



substance 


y reaction(s) Reference 


CioHg 

P1VN/PMMA film 
TMMC 


0.52 - 0.59 •• -> < 

0.47(3) •• -> < 
0.48(4) -> 


• o ! ! 

• o [IZJ 



with d < d*, as might occur for example in nanodevices, fluctuation effects 
become dominant. 

The importance of fluctuations in low-dimensional reaction-diffusion pro- 
cesses has also been confirmed experimentally. An effectively one-dimensional 
setting can be achieved by studying the kinetics of excitons (localized elec- 
tronic excitations) along polymer chains (other examples are reviewed in 
[Tot IH]). For details, consult the reviews by Kroon and Sprik and Kopelman 
and Lin in [1] . The only purpose of the polymer chains is to provide a carrier 
for the excitons. Schematically, single excitons may hop from one monomer 
to the next (thus modelling a diffusive motion) while a reaction occurs if 
two excitons meet. One may have one or both of the reactions •• — > oo or 
•• —> «o, see table ^ We shall show in section 6 that for any branching 
ratio B = r(«« — > oo)/r(«« — > »o) the long-time behaviour is still described 
by the model with a renormalized rate. For late times, one expects the 
mean exciton density to fall off as a power law n(t) ~ t~ y . The excitons are 
unstable, with lifetimes of the order < 10~ 3 s. Their decay produces light 
of a characteristic frequency whose intensity can be used to measure n(t) 
while light with a different frequency is emitted if excitons decay through a 
pair reaction. This allows to measure d t n(t) as well through time-resolved 
experiments down to the picosecond scale. Table Ogives some results for the 
exponent y in different materials (the branching ratio B < 10% in the first 
two lines of table Clearly y ~ 1/2 as expected from © and far from unity. 
This is strong evidence in favour of strong fluctuation effects in these systems 
and against their description through a react ion- diffusion equation (j2J). 

Another aspect becomes apparent if we now briefly consider the triple 
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annihilation process • • • — > o o o combined with single particle diffusion. 
The reaction-diffusion equation reads 



\n(r, t) = DV 2 n(r, t) - \n{r, t) 3 



(8) 



Following the same lines as above, but using now the Holder inequality, the 
differential inequality dtn(t) < —Xn(t) 3 leads to the bound 



This of the same order of magnitude as the long-time behaviour expected 
from diffusive motion in one spatial dimension. Therefore, and in agreement 
with scaling arguments showing that d* — 1 18J, already for the triple an- 
nihilation process, fluctuation effects should not play a major role in any 
physically realizable dimension (d > l). 2 

In conclusion, the description of reaction-diffusion processes with pair 
reactions in low dimensions requires a truly microscopic approach beyond 
kinetic reaction-diffusion equations while these equations may well turn out 
to be adequate for multiparticle reactions. For that reason, we shall in the 
following consider the master equation formulation of reaction-diffusion pro- 
cesses with pair-reaction terms only. 

2 Quantum Hamiltonian formulation 

We now review the reformulation of a non-equilibrium stochastic system 
defined by some master equation in terms of the spectral properties of an 
associated quantum Hamiltonian H and which goes back at least to the 
classic paper by Glauber [20] . To be specific, we consider in this section 
only systems defined on a chain with L sites and two allowed states per 
site. We represent the states of the system in terms of spin configurations 
u = {a i, (72, ... , <7l} where a — +1 corresponds to an empty site and a = — 1 
corresponds to a site occupied by a single particle. We are interested in 

2 This does not imply however, that models containing both binary and multisite reac- 
tion terms could not have a non-trivial behaviour. For example, the phase structure of 
the pair contact process (•• — » oo, • • o — > • • •) with single particle diffusion (»o «-> o») 
is presently controversial and under very active study |19j . 



n(t) < 




1 



r l/2 



(9) 
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the probability distribution function P(er; t) of the configurations cr. Our 
starting point is the master equation for the P(cr; t) 



w{t (t)P(t; t) - w(<t -> t)P(<t; t) 



(10) 



where w(t —>■ cr) are the transition rates between the configurations r and 
cr and are assumed to be given from the phenomenology. In order to rewrite 
this as a matrix problem, one introduces a state vector 

|P) = ^P(<7;t)|<r} (11) 

cr 

and eq. (fTUJl becomes 

d t \P) = -H\P) (12) 
where the matrix elements of H are given by 

(a\H\r) = -w(t -> cr) + 5 ct , t ^w(t -> u) (13) 

The operator H describes a stochastic process since all the elements of the 
columns add up to zero. Conservation of probability J2a- ^i "^) = 1 is 
equivalent to the relation 

(s\H = (14) 

where (s\ = (cr\ is a left eigenvector of H with eigenvalue 0. Correspond- 
ingly, H has at least one right eigenvector \s) = P s (cr) \cr) with eigenvalue 
0, that is H | s) = 0. Such a vector does not evolve in time and therefore cor- 
responds to a steady-state distribution of the system. Since in general H is 
not symmetric, this steady-state vector may be highly non-trivial. Note that 
all this is completely general and applies to any stochastic process defined 
by a master equation. With a view on the processes to be studied below one 
calls H a quantum Hamiltonian and this formulation of the master equation 
the quantum Hamiltonian formalism (see CH] for recent reviews) . The rea- 
son for this choice of language is the fact that for the processes studied below 
(and, in fact, many other processes as well) H is the Hamiltonian of some 
quantum system such as the Heisenberg XXZ Hamiltonian. The steady-state 
\s) of a stochastic system corresponds in this mapping to the ground state of 
the quantum system. The probabilistic interpretation is guaranteed by the 
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following classical result. 

Theorem 1 (Hyver, Keizer, Schnackenberg [21]): For a quantum Hamilto- 
nian H which satisfies the master equation ( fTIJ) and has (s\ = J2a- as 
a left eigenstate such that (s\H = 0, the following holds, (i) There is a 
stationary state 

| S ) = ^P s (o-)k> (15) 

such that H\s) = 0. (ii) Consider the eigenvalue problem H\n) = E n \n), 
with n = 0, 1, 2 . . .. Then 

Re E n > E = (16) 

(Hi) Let \Pq) = 1^(0)) be the initial state such that the weights of the indi- 
vidual configurations satisfy < P(<r; 0) < 1 and (s\ Po) = 1. Then for all 
times t > 0, one has 

0<P(cr;t)<l and (s\P) = l (17) 

(iv) Let H : W 1 — > M. n be a linear map such that for the elements Hij of H 
holds 

n 

#y<0 , 5Z#y = VjG{l,...,n} (18) 
i=i 

Then H is a "quantum Hamiltonian" of a Markov process described by the 
master equation KWi) . 

Time-dependent averages of an observable F are given by the matrix element 
< F > (t) = ^2F(tr)P(tr;t) = (s\ F \P) = (s\Fexp(-Ht) |P ) (19) 

a 

and we see that eq. (J 16)) means that the system indeed evolves towards the 
steady-state \s), thus -P s (cr) = P(cr; oo). 

In what follows, we shall be mainly interested in averages of particle 
numbers nj at site j and their correlators. These can be expressed in the 
quantum spin formulation in terms of the projector 
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and one-point and two-point functions of the rij are then expressed as 3 

C x {j] t) = <Uj > (t) = (s\ nj \P) , C 2 {j, £; t) =<n j n e > (t) = ( s | nfit \P) 

(21) 

Two basic situations are readily distinguished from the spectrum of H. 
If in the limit of infinite lattice size L — ► oo the lowest excited states have 
a finite gap T to the ground state energy E — 0, then the averages (|2Tj) 
will approach their steady state values exponentially fast on the time-scale 
t = 1/r. On the other hand, if there is in the L —>■ oo limit a continuous 
spectrum down to E = 0, one expects an algebraic decay of the correlators 
as the system approaches the steady state. 



3 Hecke algebra and integrability 

Before we shall write down quantum Hamiltonians for certain reaction-dif- 
fusion systems explicitly, we need some background information on Hecke 
algebras in order to make contact with integrability. The Hecke algebra 
H n (q) is spanned by n generators 

H n (q) = { ei ,e 2 ,...,e n } (22) 

which satisfy the following relations 

CiCj = Cjti ; if \i — j\ > 2 (23) 
e? = (q + <r X ) et 

where q G C is a parameter. The representations of H n (q) and the relation- 
ship to equilibrium statistical mechanics are discussed in great detail in [22] • 
We are not only interested in H n (q), but also in some quotients, denoted by 
(P, M)H n (q) (221 • Two specific examples will be of interest to us. The first 
such quotient is the celebrated Temperley-Lieb algebra (2, 0)H n (q), where in 
addition to eq. (J23|) the additional relations 

eie i±1 ei - e { = , ei± 1 e i e i ± 1 - e i±1 = (24) 

3 We stress that the structure of these matrix elements is quite distinct from expectation 
values (0| F |0) in ordinary quantum mechanics. 
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hold. On the other hand, the quotient (1, l)H n (q) is defined through the 
condition [23] 



(e;e i+2 ) e i+1 (q + q 1 - a) (q + q 1 - e i+2 ) =0 ; i = 1, 2, . . . (25) 

For the definition of more general quotients, we refer to [23J. 

Consider the N x N matrices E ab where a,b — 0, 1, . . . , N — 1. The only 
non-vanishing element of E ab is the one on the a th line and the 6 th column 
and this element is equal to unity. Define further 

E? = 1 ® ■ • ■ <g> 1 ® £ ab ® 1 ® • ■ • ® 1 (26) 

where _E ab occurs on position i and i = 1, . . . , L runs over the sites of a chain. 
Explicit realizations of the Hecke algebra or one of its quotients may be found 
in the Perk-Shultz models [25; whose Hamiltonian is of the form 

H (P,M) = ef > M) (27) 
i=i 

The importance of this observation becomes clear from the following 

Theorem 2 (Jones [2H]) : If H — Y2i=i e i' where the e,- L are the generators 
of the Hecke algebra Hi-i(q), H is integrable through the Baxterization pro- 
cedure. 

The Baxterization procedure allows to define, starting from H, in a system- 
atic way Boltzmann weights which satisfy the Yang-Baxter equation. We 
shall illustrate this in the example of the seven- vertex model in section 5. 
In many practical applications, the following result is useful. 

Theorem 3 (Martin, Rittenberg j22]) : If H — ^^T^ e^, where the are 
the generators of the Hecke algebra quotient (P, M)Hi_i{q) and furthermore 
H' = J2i=i f> where the fa are different generators of the same quotient 
(P, M)Hi-i{q) , then H and H' have the same eigenvalues, up to degenera- 
cies. 

We finish this section by writing explicit examples for the quotients (2, 0) 
and (1,1) in the case N = 2 j2E]- Then the matrices E ab can be expressed 
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through Pauli matrices 

-(;;)■- (? o).- w 

and we define a^' y,z by analogy with ()26|). Set 



e (2,0) _ _ I 



ff i + °i + — 2 ' < i+1 ~ f 2 ' 1 ~ 1+1 ' 



2 

Therefore, the Hamiltonian [2*7} EH] 



(29) 



L-l 



^(2,0) = ^ e (2,0) (3Q) 
i=l 

is integrable. In addition, it satisfies a quantum group invariance, since 

[H^°\S Z ] = [H^,^] =0 (31) 
where, recalling also that a ± = \{o x ± icr y ) 

L L 

2 



i=l i=l 

S t = «P E «*) »f «P f-^E <*) P2) 

V £=0 / V £=i+l / 

which in turn obey the commutation relations of U q (su(2)), namely 

2S Z —2S Z 

\S Z ,S ± ]=±S ± , \S + ,S-}=^- ~ 9 : (33) 

L J q — q 

On the other hand, set 

= ~~ [of^x + "Mf-i + ?«? + " <? " (T 1 ] (34) 

Here, the associated integrable Hamiltonian [2| H^ 1 ' 1 ^ = ^2f=i ef~ in ad- 
dition is invariant under the supersymmetric quantum group U q (su(l\l)), 

S1I1CG 

[#(1.1), 5"] = [#( 1 ' 1 ),T ± ] =0 (35) 
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where 



T ± =g d-D/2^ g i-l exp 



L 




± 



(36) 



and 



{T + ,T-} 



q L — q 



-L 



(37) 



-i 



4 Single-species models 

We are now ready to study explicit examples of stochastic quantum Hamilto- 
nians. The classical example merely considers particles of a single species (•) 
which may hop randomly onto an empty nearest-neighbour site (o), thereby 
modelling the reversible reaction »o <-> o» with rate D. This process is often 
called the symmetric exclusion process. The quantum Hamiltonian reads 



and coincides with the (ferromagnetic) XXX Heisenberg quantum chain [3U| . 

Certainly, one may now use the Bethe ansatz solution of Hxxx to rederive 
known results on simple diffusion. The recent interest in this setup comes 
from the insight that the integrability of the associated quantum chains al- 
low to make contact with the pre-established algebraic techniques for the 
treatment of these [23 EI] • Independently, integrability was also observed 
to occur in the transfer matrices for discrete-time dynamics [3*^1 13*3*] . The 
enormous possibilities for non-trivial applications then triggered an ongoing 
wave of activity, see e.g. |2J IH UHl I2S1 IMl ESI EHl E3 EHl ESI and references 



Following [213 EH]; we now gi ve more examples of integrable quantum 
Hamiltonians of stochastic systems, restricting ourselves for simplicity to 
a single species of particles and to binary reactions only (see section 1). 
The reaction rates are defined in table [21 using the convention of various 
authors, but unfortunately there is no standard notation. While I prefer 
a light notation (slightly modified from [40 ) and shall use it here, 4 other 

The letter v is inspired by naissance (French for birth) and a comes from Schopfung 
(German for creation). The letter [3 might have come from branching. 




(38) 



therein. 
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diffusion to the left 


o« - 


-> mo 


Dl 


«32 


W32 


Wl, i(l,0) 


1 'Ul 

1 10 


diffusion to the right 


• o - 


-> om 


Dr 


a 2 3 




^i,i(0, 1) 


-nl Q 
1 1 -LU 

1 01 


* l * 1 j * 

pair annihilation 


- 


-> oo 


2a 


0-14 


W M 


u>i,i(0, 0) 


pll 
1 00 


coagulation to the right 


- 


-> om 


lR 




W 24 


^i,o(0, 1) 


"nil 

r oJ 


coagulation to the left 


- 


-> mo 


1L 


G-34 


W34 


w ,i(l,OJ 


Til 

r io 


death at the left 


• o - 


-> oo 


h 


ai3 


™13 


wi, (0,0) 


plO 
1 00 


death at the right 


o« - 


-» oo 


Sr 


ai2 


W12 


w o ,i(0, 0) 


pOl 
1 00 


decoagulation to the left 


o« - 


- 


Pl 


042 


W 42 


^i,o(l, 1) 


pOl 
1 11 


decoagulation to the right 


• o - 


- 


Pr 


043 


W43 


^o,i(l, 1) 


plO 
1 11 


birth at the right 


oo - 


-> o« 


vr 


<*21 


W 2 1 


wo,i(0, 1) 


pOO 
1 01 


birth at the left 


oo - 


-> mo 


vl 


031 


W3I 


^i,o(l,0) 


pOO 
1 10 


pair creation 


oo - 


-> 


2a 


041 


W41 


™i,i(l, 1) 


pOO 

1 11 


Rates defined after reference 


40, 


35J 


[lOj 


[31] 


[41j 



Table 2: Two-sites reaction-diffusion processes of a single species and their 
rates as denoted by various authors. 



authors often opt for a systematic, though heavier notation with several 
indices. 

For the time being and for purposes of illustration let us consider be- 
sides diffusion only those reactions which irreversibly reduce the number of 
particles (that is, /3l,r = vl,r = cr = 0). Define 

(39) 

A = ^(q + q- 1 )(l + 6- 1 )-a/D , h = ± (2a/D + 7 (q + q~ l ) ) . (40) 

Note that the ratio of the left and right rates is taken to be the same for 
diffusion, coagulation and death processes. We first consider an open chain 
with L sites. Then the quantum Hamiltonian becomes 

H = D ( H xxz (h, A, 5)+H a + H 1 + H 5 ) (41) 

where H xxz (h, A, 5) is the standard XXZ quantum chain, including bulk and 
boundary magnetic fields 

1 

H xxz (h,A,6) = —- ^2 [ a j a j+i + a j a ]+i + A (crjcrf+i — l) (42) 

j=l 

12 




+ h fa + a* +1 - 2) - -(1 -5){q- g- 1 ) [a] - o* j+x ) 

which contains the diagonal and diffusive matrix elements while the particle 
annihilation terms are contained in 

L-1 

H a = -2aJ2r 2j ~^M+i 

L-1 

H, = -7Sg" J '(nA + ?"Vtnj + i) (43) 
3=1 

L-l 

H s = -5^q- 3 {q- 2 {l-fi 3 )aj +1 + qa+{l-fi j+l )) 

and <j ± = \{o~ x ± io~ y ) are the one-particle annihilation/creation operators. 
For a physical understanding of this we consider two special cases. 

1. Consider pure asymmetric diffusion, that is a = 7 = 5 = 0, also 
referred to as the asymmetric exclusion process. Then H = DH^ 2 ' ^ as 
given by eqs. ()29J3()j) . We thus have a very clear physical interpretation 
of the quantum-group parameter q = Dr/ Dl [3*T] . Besides simple 
biased diffusion, this model is related e.g. to the ID Kardar-Parisi- 
Zhang equation or to the noisy Burgers equation [12] ■ The quantum 
group may be used for the calculation of correlation functions [43 . 

2. In addition to diffusion, add annihilation such that 2a = Dl + Dr, 
that is A = and keep 7 = 5 = 0. Then H = D(Hq + Hi), where the 
hermitian part H = H^' 1 ^ is given by eq. (jMj) and this part alone is 
therefore supersymmetric. On the other hand, H = D X^i 1 where 
fi G (1, l)H L _i(q), but the fi are no longer symmetric. This was the 
first example of a non-symmetric realization of a Hecke algebra [HI]. We 
remark that besides the already established integrability, this system 
is also soluble through free-fermion techniques. 

Proposition 1 |3*T] : The spectrum of H in eq. is independent of the 
particle-reaction terms contained in eq. \4'J^ , that is 

spec(#) = spec(M XX z(/>, A, 5)) (44) 
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To see this, recall that the XXZ Hamiltonian conserves the number of parti- 
cles while the reaction terms irreversibly decrease the total particle number. 
Thus, H can be written in a block diagonal form 



H 



( Mo x s 
Mi 



\ 



X a 
M 2 



x n 



x. 



7,(5 



/ 



(45) 



where M n refers to the n-particle states and X are the reaction matrix ele- 
ments. Because of the identity 



det 



A X 
B 



det ^ldeti3 



(46) 



it follows that the elements of (|43|) do not enter into the characteristic poly- 
nomial of H. q.e.d. 

Therefore, the phase diagram for the full Hamiltonian H can be read off 
from the well-known spectrum of Hxxz(h,A,5) [44]. For our purposes, we 
need the following |31j . From (|40j) . only the portion of the phase diagram 
where h + A > 1 is important for us. First, he spectrum always has a 
finite gap when h + A > 1, which is realized whenever 5 ^ or q ^ 1. 
Then the ground state of Hxxz is a trivial ferromagnetic frozen state which 
corresponds to the empty state o o . . . o o. The energy gap V = E\ — Eq 
is finite. Second, the spectrum is gapless for A + h = 1, where the system 
undergoes a Pokrovsky-Talapov transition. This situation occurs for 5 = 
and q — 1. We have thus identified the cases where the model approaches 
the steady state exponentially (non- vanishing gap) or algebraically (gapless). 

At this point, it is of interest to discuss the role of the boundary conditions 
and we consider now a periodic chain, for simplicity just for the asymmetric 
exclusion process (that is a = 7 = 5 = 0). We stress that if q ^ 1, H pcr 
cannot be read from (J41|) by simply taking periodic boundary conditions. 
Rather, we have 



per 



q + q 



i=l 



q + q' 



(47) 



together with the periodic boundary conditions cr^ +1 = and cr£ +1 = af. 
The hopping terms can be brought back to the familiar XXZ form of eq. (|41j) 
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through a similarity transformation H' peT = UH pcr U 1 with the matrix 14 
exp (ng J2e=i ^ a t) w ^h q = e 2n9 such that jlSj 



#1 



1 L 
-"2(r/- fr -)f^ 



which looks the same as DH^ 2,0 \ but we now have the non-periodic boundary 
conditions <t^ +1 = q TL af, cr£ +1 — erf- As a consequence, spec(-ff per ) has 
no gap even for q ^ 1. While for a finite number r of particles and long 
chains L — > oo, it is easy to see that Re -E per ~ L" 2 |35], for finite densities 
n = r/L, an elaborate Bethe ansatz calculation shows that Re E pev ~ L~ 3//2 
^21- Therefore, and quite in distinction with equilibrium systems, a change in 
the boundary conditions may well induce a phase transition in the long-time 
behaviour (observed first in driven diffusive systems [46J). What happens is 
easily understood in this particular example. For an open chain, the particles 
get stuck at one end of the chain and a non-trivial position-dependent steady- 
state density-profile p s (i) builds up. The time needed for this should be of the 
order of the time the particles need to move from one end to the other which 
is finite for q ^ 1. On the other hand, for a periodic chain the particles keep 
chasing other forever and a steady-state particle current will be observed. By 
going to a reference frame co- moving with the mean velocity of that current, 
one is back to the case q = 1 of unbiased diffusion. 

The energy gaps T can now be found by concentrating on the spectrum- 
generating part Hxxz- For h+A > 1, the gaps are finite and are easily found. 
We concentrate here on the case of unbiased diffusion, when h + A = 1 and 
we are at the Pokrovsky-Talapov transition. The low-lying energy gaps are 
given by the following 

Proposition 2: On the Pokrovsky-Talapov line h + A = 1, 8 — 0, and for 
L large, the low-lying eigenvalues of H eq. are for periodic boundary 
conditions 

£ (P0.) = D f?A 2 { J2 + + J2) _ **D A £ (Jj _ ^ + Q (i _ 4) 

(49) 

where the Ij are pairwise distinct integers (half-integers) whenr is odd (even). 
For an open chain 

E " me) = d {l) 2 ^ 1 + --- + ® • ( x " 2( 7 1} i A a ) + ^ (50) 
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where the Ij are pairwise distinct non-negative integers. The integer r = 
0, 1, 2, . . . gives the number of particles in the sectors of H X xz- 
The finite-size amplitudes lim^oo L 2 E r are independent of A, for either 
periodic or open chains. This proves an old conjecture [31] based on numerical 
calculations. For A = 0, the well-known free-fermion solution is reproduced. 
To leading order in eigenvalues with the same ■ -+1? are degenerate. 
We observe that for periodic boundary conditions, this degeneracy is already 
broken by the first correction in 1/L, while for free boundary conditions, the 
leading correction keeps that symmetry. 

Eqs. f|49!50p are easily found from the Bethe ansatz. We first consider 
periodic boundary conditions. The XXZ chain may broken into sectors con- 
taining only the states with r particles. Performing the Bethe ansatz as usual 
[4*Tj . one has for the energies 

E r = 2D (r — cos k\ — . . . — cos k r ) (51) 

where the quasimomenta k±, . . . , k r are solutions of the Bethe ansatz equa- 
tions 

r 

Lh j = 2*I j -^Q(k j ,k l ) ; j = l,...,r (52) 
i=\ 

where the Ij are pairwise distinct integers (half-integers) when r is odd (even) 
and 

r^M ,a „ Asin((A; - k')/2) 

Q(k, k ) = 2 arctan — , w " . — ^777 r--- 53 

V ' 1 cos((k + k')/2)-Acos((k-k')/2) K ' 

We are interested in the leading finite-size corrections when L —>■ 00 with r 
fixed. The ansatz 

2tT Q/j ,_ JN 

kj = — J; + ^ + . . . (54) 

gives Q(h,k 2 ) ~ (Ji-J^ + O^- 2 ) and a, = g£ ELi(^~ h)- Th en 

(J4TJJ) follows. Second, for free boundary conditions, the Bethe ansatz [IS] 
reproduces eq. (foTj) for the energies, while the Bethe ansatz equations for the 
quasimomenta now take the form 

Lk 3 = nlj - 1 £ [Q(kj, ki) - Q{-k j} ke)] (55) 

for j — 1, . . . , r and where the Ij are pairwise distinct non- negative integers. 
Eq. (J53J1 leads to aj = —j^(r — l)irlj and we arrive at (j5UJ). q.e.d. 
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5 The seven- vertex model 



Having seen for some single-species reaction-diffusion processes how the re- 
lationship with the Bethe-ansatz solution XXZ chain could be used to infer 
certain physical properties, we present in this section an example of the Bax- 
terization procedure That procedure permits to associate to a stochastic 
quantum Hamiltonian H related to a Hecke algebra the Boltzmann weights 
of a corresponding two-dimensional vertex model and thus prove its integra- 
bility. In principle, the model is then solved through the Bethe ansatz. For 
the six- and eight-vertex models the completeness of the Bethe ansatz has 
recently been proven 

Following (HI], we consider the pair-annihilation model eq. (J41|) already 
defined in section 4 with 7 = 5 = and furthermore, we take 2a = Dr + Dl, 
thus A = 0. We call Q = 2a/ D. After having performed the canonical trans- 
formation Ef -> (-l) a - b Ef , only at even sites i, the quantum Hamiltonian 
takes in the basis given by eq. (|28|) the simple form 



L-l 



H 



(56) 



i=l 



where 



l x ® ...<g) ® 



/0 







r 1 



o 








Lj+2 







(57) 



q + q- L j 



and lj are 2x2 unit matrices attached to the site i. 

While we could certainly solve this particular model through a Jordan- 
Wigner transformation followed by a canonical transformation, see e.g. |4T7] . 
we are interested in generic approaches of a more general validity than in 
those cases reducible to a free-fermion description. 

We have already seen in section 4 that the satisfy the Hecke algebra 
(|23|) . We now construct a two-dimensional vertex model corresponding to H 
having a row-to-row transfer matrix T(6) depending on the spectral para- 
meter 9. That transfer matrices will satisfy the Yang-Baxter equations [5I| 
which imply the commutation relations [T(8),T(6')] = if 9 ^ 6'. The 
construction is based on the matrix Ri(9), % = 1, 2, . . . , L — 1 which depends 
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on the spectral parameter 6. The Baxterization procedure for Hecke algebras 
[25] gives 

- ,„ s sinh6> sinhf?? — 0) „ . , 

R ^ d ) = —T~ e * + H " : 9 = e" 58 

smh r; sinh 77 

which for our model (|5fi|57|) leads to 



Ri(0) = li 



1; 



i-1 



72. 



M+i 



1; 



i+2 



(59) 



with the non-vanishing elements of the matrix 72j j + i 



sinhr/ 



/ sinh(r7 



V 



e 61 sinh 77 sinh 
e~^sinh0 e 6 * sinh 77 



ft sinh \ 



sinh(r/ + 9) J 



(60) 



The relations eq. (jzHj) imply that these matrices satisfy the spectral parameter- 
dependent braid-group relations 



Ri{9)R 



i±i\ 



[Rilfi)M 



; |i-j'|>2 



(61) 



which are equivalent to the Yang-Baxter equations. 

In a 2D vertex model with vertex configurations labelled by (k, i, m, n) 
the Boltzmann weights Sj^ are obtained from 



Ri(0) 



Qk,n -t 



4-1 



E 



mk 



E 



nl 



4+2 



(62) 



This implies that the vertex model associated to eq. (J51)Jl is a seven-vertex 
model. In a vertex model, arrows are attached to the bonds of a square 
lattice jHU]. In the stochastic model, we associate a particle (•) with an arrow 
pointing up/right and no particle (o) with an arrow pointing down/left. In 
figure El we list together the chemical reactions, the vertex configurations and 
their Boltzmann weights. The vertices usually labelled 1 to 4 correspond to 
no reaction and are not shown (see [31j). Vertices 5 and 6 correspond to 
diffusion to the right and to the left and vertex 7 to pair-annihilation. In the 
leftmost column of figure El the state of the particles before the reaction is 
given as the lower pair of symbols while the state after the reaction is given by 
the upper pair of symbols. The middle column gives the corresponding vertex 
configuration and the right column the Boltzmann weight. The Hamiltonian 
eq. may be recovered from H = — ■^ha.T(9)\ e . 
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Figure 2: Diffusion and pair- annihilation of particles in the seven- vertex 
model and their Boltzmann weights, for the vertices 5 to 7. 



5. 




e v sinh 9 



6. 




e v sinh 9 



7. 




Q sinh 9 



6 Further applications 

We have studied in some detail the pair-annihilation process and its integra- 
bility. Still, extracting explicitly information about the long-time behaviour 
(or the steady-state in more complicated models) is not yet trivial. In this 
section, we briefly review some approaches which may be useful. 



A) Spectral integrability. We have already seen that in certain cases, 
the quantum Hamiltonian H — H + Hi such that spec(if) = spec(ifo), 
independently of the precise form of Hi. It may happen that although H is 
integrable, H is not. Such a model is said to be spectrally integrable. If only 
binary interactions are present, it is convenient to express H in terms of a 
two-site matrix H i i+ i acting on the sites % and i + 1 

L 

H = li ® • • • ® lj-i <8> H jJ+ i ® l j+2 ■ ■ ■ <g) 1 L (63) 

3=1 
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s 


Ci(t) C 2:n (t) 




< a + 7 
> a + 7 


r l/2 £-3/2 

exp(-25t) r 1 / 2 exp(-45t) 
exp(-2<ft) exp [-4<ft + 2(A + 1/A - 2)t] 



Table 3: Generic long-time behaviour of the one-point and two-point cor- 
relators C\(t) and C2- n {t) for n finite, in the system eq. (|64|) and with a 
translation-invariant initial state and a finite initial particle-density. 



For the model (}4*Tj) with left-right symmetry, that is q — 1, we have 



[0 


-5 


-5 


-2a \ 





1 + 5 


-1 


-7 





-1 


1 + 5 


-7 










2(a + 7 ) / 



(64) 



One can always rescale time such that D — 1. Then the parameters of the 
XXZ chain H = H X xz become A = 1 + 5 — 7 — a and h = a + 7. 

In eqs. (J63I64|) is only integrable for 5 = 0. On the other hand, the 
special case a = 0,7 = 5 simply corresponds to the radioactive decay of 
diffusively moving particles. While for 5 = 0, one is back to the critical 
Pokrovsky-Talapov line h + A = 1, the associated quantum spin chain has a 
frozen ground state for 5 = 0. The one-particle and two-particle correlators 
= X^-=iCiO';*) and C 2;n (t) = ^ =1 C(j,j +n;t) with n fixed, see 
eq. (}2*Tj) . only imply the sectors with r = 1 and r = 2 particles of i^xxz, 
respectively. Their long-time behaviour is easily worked out from the results 
of section 4 and is collected in table EH 4Qj. Of course, we implicitly assume 
that the corresponding amplitudes do not accidentally vanish. At first sight, 
one might have expected a simple exponential factor e~ 2k6t for the /c-point 
correlator Ck and we already observe that eventual algebraic prefactors are 
not readily predicted from the spectrum of H alone. The more complicated 
form of the relaxation time for 5 > a + 7 comes from a bound state in the 
two-particle sector of Hxxz, with energy 45 + 4 — 2 A — 2/ A, see (HDl EH]- 
More general initial conditions are discussed in |4Uj . 

B) Similarity transformations. In trying to extract explicit information 
on certain reaction-diffusion systems, the integrability of the quantum Hamil- 
tonian H plays a central role. Since it is difficult to realize the constraints 
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(|18|) of stochasticity and integrability at the same time, it is of interest to 
see whether there exist systematic transformations of an integrable quantum 
Hamiltonian H towards a new stochastic Hamiltonian H. 
Specifically, we shall consider the transformation 

L 

H = BHB~ 1 , B = (g)B J (65) 

i=i 

where Bj = 1®...®1®5®1<S>...<8>1 is the transformation matrix B 
acting on the site j. Then the systems described by H and H are said to 
be similar to each other. An interesting alternative, which has not yet been 
systematically studied, is to consider an enantiodromy transformation JU| 

H = BH T B~ l (66) 

where H T is the transpose of H. 

From now on, we shall focus on translationally invariant systems and 
consider periodic boundary conditions. The effect of the transformation B 
on H is completely given by its effect on the two-particle Hamiltonian Hjj + % 
in (|64|). A stochastic similarity transformation arises if both H and H rep- 
resent stochastic systems. For a simple example, consider the symmetric 
annihilation-coagulation process eqs. (|63I64|) with 5 = 0. If C(t\a, 7) = C\(t) 
is the spatially averaged particle density with the rates a and 7, respectively, 
a stochastic similarity transformation shows that EH1 EI] 

C x (t\a, 7) = ir~~~r~~ Ci(*|0, a + 7) (67) 
la + 7 

Similar results hold for any /c-point correlator Ck(t). So far, explicit methods 
to find the time-dependent correlators are only available for either the pure 
coagulation model a = through empty-interval methods (see below) or 
the pure annihilation model 7 = through free-fermion techniques, see |55| 
E3 EI] and references therein. Eq. (jfi7j) allows to reduce any symmetric 
annihilation-coagulation process to pure coagulation, for any initial density 
Ci(0|0, 7). This also explains the experimental results in tabled The known 
stochastic similarity transformations of the form (|65J) leave the parameters 
A and h of the XXZ chain invariant, but the results of Proposition 2 suggest 
that a stochastic similarity transformation between systems with different 
values of A might exist. See (HS] for the extensions to 5 > and q ^ 1. 
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model 


reactions 


conditions 


A 


• • <— > oo 






• o <- 


o« 


2(a + a) = D L + D R 


B 


• • — > oo 


— > o«, 


• o 


• o <- 


o« 


2a + 7l + 7h = D L + D R 


c 


oo — > 


oo — > o«, 


• o 


• o <- 


o» 


2a + u L + u R = D L + D R 


D 


— > mo o» 


• o, o« — > 




• o <- 


-± o« 


lL = D L J R = D R 


E 


OO — !> wo, o» 


• o, o« — > 


oo 


• o <- 


o« 


V L = D L Vr = Dr 


F 


• o, o» — > oo 


• o, o» — > 








5r/5l = Pr/Pl 


G 


• o, o« <- 


-» ••. oo 








f S R = Pr 7ij = u L 

\ $L = Pl lL = Vr 


H 


• — s> o 


o — S> • 








f h = $R = IL = 1R 
\ Pl = Pr = V L = Vr 



Table 4: Single-species processes with space-independent reaction rates and 
which are similar via (|65|) to a free-fermion model. The reaction rates are 
defined in tabled 



Eq. (joTj) allows to recover the long-time behaviour C(t\a, j) ~ t -1 / 2 from 
a simple heuristic argument. For pure coagulation, one particle always re- 
mains, thus C(oo|0,7) = 1/L in the steady state. Therefore the steady-state 
density C(oo|a, 7) ~ L" 1 . On the other hand, from the spectrum of H, the 
leading relaxation time r = T -1 ~ L 2 ~ £ 2 , where £ is identified as the 
characteristic spatial length scale. Therefore C(oo|a,7) ~ ~ r -1 / 2 . The 
asserted time-dependent behaviour therefore might have been anticipated on 
dimensional grounds. 

C) Free fermions. For the pure annihilation model with 2a = Dr + Dl, 
one has A = 0. In this case, H may be diagonalized through a Jordan- 
Wigner transformation followed by a canonical transformation jJS]. In order 
for this to work, H may only contain pairs of particle creation and annihi- 
lation operators. For space-independent reaction rates, the complete list of 
react ion- diffusion process whose quantum Hamiltonian H is similar through 
()65|) to a free-fermion Hamiltonian H is as follows jTOJ HTJJ E3] and shown 
in table 0] Since the transformation (|65|) is spatially local, these correspon- 
dences actually hold in any dimension, but free-fermion methods are only 
available in ID. 

Of these models, only models A (diffusive pair-annihilation and creation, 
solved exactly in ID in G (kinetic Ising model with Glauber dynam- 

ics) and H (free decay and creation of particles) are reversible and have an 



22 



equilibrium steady-state. Their quantum Hamiltonian is therefore similar to 
a symmetric matrix. The similarity of the kinetic Ising model with Glauber 
dynamics (model G) to a free-fermion model was obtained long ago through 
a duality transformation jSH] and more recently as a true similarity transfor- 
mation EI]. This suggests to study a more general type of relationship, 
based on domain-wall dualities, see (TU1 EZj for details. Models C and E 
are obtained by a particle-hole permutation • <-> o from models B and D, 
respectively. Model B is the biased annihilation-coagulation process, while 
model D is the diffusive coagulation process with arbitrary decoagulation. 
Finally, model F is the doubly biased voter model (in space and in the pref- 
erence between • and o) and in j^Hj some correlators are found from the 
free-fermion form. The physical behaviour of all these models can be treated 
in a single calculation. For example, the mean particle-density depends on a 
single parameter h such that |3U] 

f r 1 / 2 • i£h = l 

w v ' \ t • j / i exp(-t/r) ; if < h < 1 v ' 

where Ci(oo) is the steady-state density and r = 1/(4 — Ah) is the rela- 
xation time (see JU] and references therein for more information on solved 
free-fermion models). 

This kind of analysis was generalized to find those reaction-diffusion sys- 
tems which are similar, via a transformation of the type (|65j) . to the XXZ 
chain j^UJ. While the full result is too complex to be re-stated here, an 
interesting special case is given by the conditions 

-f R + (3 L + 2a + D L = u L + 5 L + 2a + D R 

lL + f3 R + 2a + D R = u R + 6 R + 2a + D L (69) 

In this case the (usually infinite) hierarchy of equations of motion for the k- 
point particle-density correlators Chit) = < riiit) . . . rikit) > closes naturally, 
such that Ck(t) only depends on the Ce(t) with £ < k. In principle the equa- 
tions of motions for the Ck can then be solved iteratively [53] . 

D) Partial integrability. The previous sections have shown that con- 
structing integrable stochastic systems which go beyond mere free diffusion 
is a non-trivial exercice. One might wonder whether the condition of full in- 
tegrability is not too strong. After all, from a practical point of view it would 
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be enough to identify a set {Qi, . . . , Qm} of observables such that these sat- 
isfy a closed set of equations, say Qi = fi(Q%, ■ ■ ■ ,Qm), with i = 1, . . . , M. 
Such a partial integrability may be enough for many practical needs. Indeed, 
such an approach is available through the empty-interval method [12 Ej- Con- 
sider a periodic chain with L sites and lattice spacing a. Let I n (t) be the 
probability that at time t, n consecutive sites are empty Then the mean 
particle-density is [12 E2] 

C 1 (t) = (l-I 1 (t))/a (70) 

In order to illustrate the method, we consider the left-right symmetric pure 
coagulation model and also take the free-fermion condition 7 = D of model 
D in table but we now add a three-site production reaction • o • 
with rate 2DX [SH]. The equations of motion for the I n (t) read 

h(t) = 2D(I (t)-2I 1 (t) + I 2 (t))-2DX(I l (t)-2I 2 (t) + h(t)) 
i n (t) = 2D (I n ^(t) - 2I n (t) + I n+1 (t)) ; 2<n<L-l (71) 

together with the boundary conditions Io(t) = 1 and I hit) = (assuming 
that there is at least one particle in the system). The solution of these equa- 
tions is straightforward. For example, one obtains for the leading relaxation 
time r" 1 = T — 2Dtt 2 L~ 2 + 0(L~ 4 ), in agreement with the results of sec- 
tion 4. The effect of the production term is only transient, as illustrated in 
figure El for the mean density C\(t). For A = 0, C\(t) ~ l/Vt is of course ex- 
pected from eqs. (J6I67|) . While the free-fermion condition 7 = D is essential 
for the method to work, we also see that the presence of the production term 
poses no problem at all for the closure of the equations of motion (|7Tj) . 5 

Accepting the free-fermion condition j LjR = D L Rl one can extend the 
treatment to the more general model D of table 0] and may even extend 
this further to include the processes 00 — > and 00 — > o«, «o with rates 
2cr, v R , v L , respectively [3H2]- Let us call this system model D ? which de- 
pends on the seven parameters D L R , Pl jR , vl,r, o~. In a remarkable paper [4T] . 
the idea of the empty-interval method was translated into the Hamiltonian 
formalism and several new sets of observables were defined which generalize 
the variables I n (t) and lead again to closed equations of motion. It turns out 
that the spectrum of relaxation times of model D' is given by the Hamiltonian 

°This term cannot, e.g. by a similarity transformation, be turned into a term treatable 
by either free-fermion or full integrability methods. 
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Figure 3: Evolution of the mean particle density C\(t) in the symmetric 
coagulation model with the production reaction • o • — > • • • for several 
production rates A. For long times, the asymptotic behaviour is Ci(t) ~ l/y/i 
for all values of A (after 

of the Wannier-Stark ladder [UJ 

L 

H = -H + « + i + ( h + h ' n X] (72) 

n=—L 

where h and h' are constants. In this case, the couplings in H are space- 
dependent. The extension of the similarity/enantiodromy approach to this 
more general setting remains to be done. Extensions of the empty-interval 
method to interactions on more than two sites are studied in [UU] . 

While the empty-interval method as such does not work for the pair- 
annihilation process, the method has been generalized recently [HI]. We 
briefly explain the idea using the left-right symmetric pair-annihilation pro- 
cess with the free-fermion condition a = D (model A or B) as example. Let 
G n (t) be the probability that at time t, one has on n consecutive sites an 
even number of particles. The mean particle density is C\{t) = (1 — Gi(t))/a. 
Furthermore, let F n (t) (H n (t)) be the probability that a segment of n consec- 
utive sites with an even (odd) number of particles is followed by the presence 
of a particle at the (n + l) th site. From the relations 

2F n (t) = (1 - G x ) + (G n - G n+1 ) , 2 J H n (t) = (l-G 1 )-(G n -G n+1 ) (73) 
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and the boundary condition Go(*) = 1; the equations of motion 

G n (t) = 2D (F n „! - H n ^ + H n - F n ) = 2D (G n _ x (t) - 2G n (t) + G n+1 (t)) 

(74) 

follow. They can be solved by standard methods. Reaction terms parame- 
trized by a, vl,r, Pl,r (see table EJ) and even the reaction o • o — > • • • can 
be added [HI]. Correlators are studied in 

In view of the practical success of these techniques it is perhaps not com- 
pletely futile to ask whether there might a be a systematic way to identify 
these 'empty-interval' or related variables ? 

E) Multi-species models. We consider chains with N states per site. One 
of them is taken to be the empty site (o) and the other states are labelled 
A n , n — 1, . . . , N — 1. Finding integrable stochastic systems becomes more 
difficult when N increases. Several examples were found through the 
quotients (P,M)Hl-x(q) as realized through the Perk-Shultz model. They 
are collected in table |5J The following conventions apply. 

1. For the first reaction in all models and the second reaction in models 
03, (£, D (with n < m understood) the reaction to the right (left) occurs 
with rate T R (T L ). 

2. For models (£, 5 the sum r + s has to be taken modiV. If in this case 
r + s = mod N, the rate is Tl + T/j. If r + s ^ mod N as well as 
for the third reaction in models £, D the rate is r#. In model D, it is 
also assumed that in the third reaction, pairs (r, s = r ± 1) never have 
an element in common. If the products on the right are interchanged, 
(e.g. AiAi — ► Aio in model <£), the rate is Tl. 

3. In the third reaction in model £ the rates are T±, respectively such 

that r + + r_ = r L + r R . 

One defines D = vT^Tr = 1 and q = \/Vr/Tl- The Hecke algebra quo- 
tient (P, M)H L _i(q) according to the realization as a Perk-Shultz quantum 
chain eq. (|77|) [2H I2E1 is also indicated. 

From tableland Theorem 3 we see that the simple diffusion model 21 has, 
up to degeneracies, the same spectrum as the XXZ chain used in section 4 to 
described biased diffusion of a single species of particles •. In the same way, 
the spectrum of model <E is, up to degeneracies, the same as the one found 
for pair-annihilation in section 4, with 2a = D^ + Dr. 
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model 


reactions 


quotient 


21 


v4 n O 


<— > 










(2,0) 




A n O 


<— > 




A A 


<— > A m A n 




(N,0) 


€ 


A„o 


<— > 




A A 

- rl n- r '-m 


<-> A m A n 




(N-1,1) 




A n O 


<— > 


Ov4„ 


A A 


^ A m A n 




(N-2,2) 




A n O 


<— > 


oA n 


A r A s 


—> oA r+s 




(1,1) 


5 


A n o 


<— > 


oA n 


A r A s 


— > oA r+s 


A n A n > A n -i-iA n -i-i 


(2,1) 



Table 5: Some integrable reaction-diffusion processes of iV — 1 species and 
their Hecke algebra quotient j2E], see text for the definition of the rates. 



For illustration, we briefly consider model (E for N = 3. Each site may 
contain either a particle of type A or B or be empty (o). Single particles 
may diffuse to the right Ao — ► oA, Bo — > ol? with a rate T# or similarly 
to the left with rate i\. On encounter, between like particles the reactions 
AA — > oB and — > Bo occur with rates Tr and I\ respectively and 
similarly BB — > oA, Ao. Two unlike particles react — > oo with rate 
Tl + T_R- In the left-right symmetric case, the identity of the spectra of 
and the one of pair-annihilation, up to degeneracies, was checked directly [3T] • 
Furthermore, in the spirit of the empty-interval method, a closed system of 
equations of motion was found, whose solution leads to the mean particle- 
densities n A (t) ~ n B (t) ~ t~ x l 2 

In [HI], Bethe-ansatz solutions of the master equation for iV-species mod- 
els with particle-numbers conservation are studied. In particular, model 03 
with iV = 3, 5 was rediscovered. The models in [S3] are found from solutions 
of quantum Yang-Baxter equations. Further study might reveal a relation- 
ship to diffusion algebras [3U], see below. 

For periodic boundary conditions and N > 2, the diffusion bias leads after 
a similarity transformation to a generalized Dzialoshinsky-Moriya interaction 
A sufficient criterion for integrability was derived in an attempt to look 
more systematically for integrable many-species models [fJo] . 

Finally, a different generalization from section 4 is to consider integrable 
stochastic models on ladders [71], rather than chains. 

F) Diffusion algebras. For certain integrable systems, there exist algebraic 
methods which allow to find the steady-state |s) such as the celebrated matrix 
product states j21E7J- Time- dependent problems are treated in 
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Behind this seemingly technical and ad hoc method there is a new and 
general mathematical structure. We shall explain here the main idea us- 
ing reaction-diffusion systems with N states per site labelled by A n , n = 
0, 1, . . . , N — 1 (where Aq = o) moving on a periodic chain with L sites, but 
generalizations to different boundary conditions are possible. The allowed 
reactions are A n A m — > A m A n with rate g nm (in particular, model 23 from 
table 03 is a special case of this). The un-normalized steady-state probability 
distribution is [2] 

P s (ct) = P(a u . . . , a L ) = Tr iV ai V a2 ■ ■ -V aL ) (75) 

where the matrices T> a satisfy the quadratic relations 

g ap V a V p - g pa V p V a = x p V a - x a V p (76) 

where a < p and a, p G {1, . . . , iV} and g ap G IR\{0}, g pa G M. and the x a are 
complex parameters. If in addition the set A of generators T> a admits a linear 
PBW basis of ordered monomials V^V^ 2 ■ ■ ■V'l" with en > <t 2 > . . . > u n 
and k n G N, A is called a diffusion algebra 

These conditions imply certain constraints on the g ap and the x a , quite 
analogously to the Jacobi identities of a Lie algebra. These constraints can 
be fully solved and a classification of all diffusion algebras for N species is 
obtained [123 EH]- The representation theory of iV-species diffusion algebras 
is just getting started, see [70] . 

7 Outlook: local scale invariance 

We finish with a discussion on how the scale invariance of many reaction- 
diffusion systems might be turned into a dynamical symmetry. For example, 
the symmetric pair-annihilation process is on the Pokrovsky-Talapov critical 
line. One has the covariance 

< n(n, h) . . . n(r p , t p ) > = b~^ + - + ^ < n(r[, t[) . . . n(r' p , t' p ) > (77) 

of the p-point correlators under the dilatation r — > r' = br, t — ► t' = b z t 
of the space and time coordinates r, t respectively, where z is the dynamical 
exponent and scaling dimensions. In the cases at hand, z = 2. 

This is reminiscent of the situation at equilibrium critical points. In 
those systems, it is known that under fairly general conditions, the covari- 
ance of the p-point correlators under global scale transformations r —>■ br can 
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Figure 4: Snapshot of the coarsening of ordered domains in the 2D Glauber- 
Ising model, after a quench to T = 1.5 < T c from a totally disordered state 
and at times t = 25 (left) and t = 275 (right) after the quench. 

be extended to conformal transformations. In addition, in two dimensions 
the energy-momentum tensor of local conformally invariant field theories be- 
comes an analytic function T = T(%) of the complex coordinate 3 such that 
not only T($) itself, but also all powers T n ($), n = 1, 2, 3 . . ., are conserved 
[Z2HZH1- This signals the integrability of 2D conformally invariant field the- 
ories. Is it possible to generalize the space-time dilatations encountered for 
critical reaction- diffusion systems in a similar way ? 

This question has been recently addressed in the context of kinetic spin- 
systems [731 ESI ■ We have already seen above that the kinetic Ising model 
with Glauber dynamics [20] m ay be obtained through a similarity transfor- 
mation of the quantum Hamiltonian from a certain single-species reaction- 
diffusion system, see model G from table EJ We now concentrate on this sys- 
tem. In the Glauber-Ising model the transition rates in the master equation 
are chosen such that the steady-state \s) is given by the equilibrium probabil- 
ity distribution P s {cr) ~ e~ n<l(7 ^ T with the classical Ising model Hamiltonian 
Tt = — a i a j where T is the temperature. Glauber dynamics may be 
realized through the discrete-time heat-bath rule o~i(t) — > ai(t + 1) such that 

<Ti(t + 1) = ±1 with probability \ [1 ± tanh(/i;(t)/T)] (78) 

with the local field hi(t) = h + ^2jU) cr j(^)- With the choice (|78j). the master 
equation can be solved exactly in ID [20] • The time-dependent spin-spin 
correlators and their approach towards equilibrium are thus determined. 
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In contrast to equilibrium statistical mechanics, where fine-tuning the 
model parameters is needed to reach a critical point, dynamical scaling is 
often found to occur in large regions of the model's parameter space. For 
example, prepare the system initially in a disordered state and then quench 
the temperature to a final temperature T < T c below the critical temperature 
T c > 0. 6 Although the steady-state of the model is not critical, the relaxation 
towards it occurs through domain coarsening and is very slow, the typical 
length scale varying with time as L(t) ~ t 1 ^, see figure 0] Typically the 
observables depend algebraically on the time t passed since the quench, see 
[HI HI] for (a collection of) recent reviews. Here we concentrate on the two- 
time spatio-temporal response function R(t, s; r) of the time- dependent spin 
a r (t) at site r with respect to an external magnetic field h (s) applied at 
the origin at an earlier time s < t. Generically, two-time quantities such 
as R(t, s; r) depend on both times t and s and not merely on the difference 
t = t — s. This breaking of time-translation invariance is called ageing. 

For ageing systems, an extension of dynamical scaling is possible and 
allows to fix the form of the two-time response function. Specifically, it can 
be shown that for a dynamical exponent z = 2 [77J EH HH] 



R(t,s;r)- SMt)) 



5h Q (s) 



h=0 



, s l+a-X R /2 

r (-J (t-s)- 1 -°exp 



M r 2 



(79) 



Here a and \r are non-equilibrium exponents to be determined which char- 
acterize the ageing universality class [TS]. Finally, r and M. are non- 
universal constants. We first present evidence that the response function 
of the Glauber-Ising model in 2D and 3D is indeed given by (|79|). Then we 
discuss where this presumably exact result comes from. 

We first consider the autoresponse R(t,s) = R(t,s;0). While R itself is 
too noisy to be measured directly, integrated response functions are accessible 
through simulations, see [JJJ] and references therein for the details which we 
skip over here. For the example, the integrated autoresponse 

p(t,s)= [ duR(t,u)~ s - a f M (t/s) (80) 

is relatively easy to measure, whereas the scaling function /m(^) can be 
calculated explicitly from (fTHj) . In the Glauber-Ising model, the exponent 

6 In the ID Glauber-Ising model, T c — leads to certain modifications of the ageing as 
described from the point of view of local scale invariance |7ri| . 
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Figure 5: Scaling form of the integrated magnetic response in the Glauber- 
Ising model as a function of x = t/s below criticality. The symbols corre- 
spond to different waiting times s. The integrated autoresponse is shown 
(a) 2D at T = 1.5 and (b) 3D at T = 3. An example of the integrated 
spatio-temporal response in 2D at T = 1.5 and with p = 2 is shown in (c). 
The full curves are obtained from ()79jl . After [7H] . 

a — 1/z — 1/2 (see jHU] for a detailed discussion) and ~ 1.26 and 1.6 in 
2D and 3D, respectively. In figure E^b, the scaling function /at (a?), as ob- 
tained from large-scale simulations, is shown for several values of the waiting 
time s. In both two and three dimensions, a nice scaling behaviour is found 
and the form of the scaling function agrees very well with the prediction 
from eq. (fTUj). Next, the r-dependence of R(t,s;r) is tested by measur- 
ing the spatio-temporally integrated response f^du J^ 1J drr d ~ 1 R(t,u;r) ~ 
s d ' 2 ~ a p ( 2 '(t/s,p), where p is a control parameter. We stress that the scaling 
function p^ does not contain any more free non-universal parameter [7H] . 
As an example, we compare in figure |U data from 2D taken with \x = 2 
with eq. (|75|). Besides the expected scaling, the functional form of the scal- 
ing function neatly follows the prediction. We stress that the position, the 
height and the width of the maximum of p^ 2 ' in figure Efc are completely fixed. 
Similar results have been obtained for other values of p and in 3D as well. 
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This provides strong evidence that eq. (|79|) is exact, at least in this model 
[7Uj . Tests of (J79|) in different universality classes are described in [75] . 
In order to derive (J7HJ), consider the diffusion equation 

(2Md t - d r ■ d r ) <f>{t, r) = (81) 

For fixed Ai, the Schrodinger group is the maximal invariance group on the 
space of solutions of eq. (JHTj) . It is defined by the space-time transformations 
(7Z is a rotation matrix) 

. at + 3 . IZr + vt + a „ . , 

and acts projectively on the solutions (f>(t, r) |8*T| . Let scf) d be the Lie algebra 
of ()82j). Time-translations occur in sct) d and are parametrized by f3. If we 
treat the 'mass' Ai not as a constant but as another variable, the embedding 
sd) d C conf d+2 for the complexified Lie algebras follows [82J, where conf d+2 
is the Lie algebra of the conformal group in d + 2 dimensions. From the 
classification of the parabolic subalgebras of conf d+2 we obtain several new 
subalgebras, called age or alt [52]. For the ID case, we illustrate in figure |H1 
their definition through the root space of conf 3 = B 2 . These subalgebras still 
contain the generator for the dilatations t — > b 2 t, r —>■ br (which is in the 
Cartan subalgebra of conf 3 ) but do not contain time-translations anymore 
(which is in the lower left corner of figure E^bc). They are candidates for a 
dynamic symmetry algebra of ageing systems. If we assume that the two-time 
response function transforms covariantly under the action of either age or a ft, 
a set of linear differential equations for R(t, s; r) is obtained. Matching their 
solution with the expected [IE] scaling behaviour of R, we recover eq. f7H|) 
in the special case z = 2. 

The functional form of R depends on the fact that the Galilei transfor- 
mation of (|H2*|) is identical to the well-known one of a free particle. It is not 
trivial at all that the response function of an interacting field theory such 
as the Glauber-Ising model in d > 1 dimensions should be recovered from a 
dynamical symmetry of the equation of motion of a free-field theory. 

There exist infinite-dimensional Lie algebras which contain scfj d as subal- 
gebras. For example, the Schrodinger group (J%2~|) is a subgroup of the group 
defined by the transformations t t' and r — ► r' where 



if = 3{t) , r' = rJd{t) or else t' = t , r' = r - a(t) (83) 
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Figure 6: Root space of the complexified conformal Lie algebra conf 3 , in- 
dicated by the full and the open points. The double circle in the center 
denotes the Cartan subalgebra. The generators belonging to the three non- 
isomorphic parabolic subalgebras [H2] are indicated by the full points, namely 
(a) scf)!, (b) oqz 1 and (c) alii. 

and (3 and ex. are arbitrary functions. Whether this has a bearing on the 
ageing behaviour of non-equilibrium spin-systems is still open. Local scale 
transformations generalizing the Schrodinger group (J%2*J) to general values of 
the dynamical exponent z 7^ 2 exist [ZS]- It can be shown that R(t, s; r) = 
R(t, s; 0)$(r(t— s) _1 / z ), such that eq. (f?H|) holds for the autoresponse R(t, s; 0) 
if \r/2 is replaced by \r/ z and $(t>) is given as the solution of a linear dif- 
ferential equation of fractional order [73]. 

I thank M. Pleimling and J. Unterberger for their pleasant collaboration on 
local scale invariance. This work was supported by CINES Montpellier (pro- 
jet pmn2095) and the Bayerisch-Franzosisches Hochschulzentrum (BFHZ). 
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